Square lattice site percolation thresholds for complex neighbourhoods 
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In this paper we compute the square lattice random sites percolation thresholds in case when sites 
from the 4th and the 5th coordination shells are included for neighbourhood. The obtained results 
support earlier claims, that (a) the coordination number and the space dimension are insufficient 
for building universal formulae for percolation thresholds and (b) that percolation threshold may 
not decrease monotonically with lattice site coordination number. 

PACS numbers: 64.60 
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I. INTRODUCTION 

The percolation [EH El 

is one of central topic of statis- 
tical physics with many applications from theory of the 
fluid flow through the porous media 0| to sociophysics Q 
via immunology @, 0] , fires spread Q , and many others 

Probably in all these fields the percolation threshold p c 
plays the crucial role: When the sites occupation proba- 
bility p is larger than this critical value p c the liquids may 
flow from one side of porous rock to another, epidemic 
or disease may spread over all individuals of a popula- 
tion, and fire may destroy all trees in a model forest, etc. 
Thus evaluation of the percolation threshold for many 
networks 0, [j| and lattices ES EE E1E1 El and even 
surfaces [l5| still appeals to scientists [l6| since the begin- 
ning of percolation theory [l| . As exact values of percola- 
tion threshold for the site percolation are known only for 
several lattices — for example: ID chain of sites (p c = 1) 
[H, the triangular lattice (p c — 1/2) [3j, Kagome lattice 
(p c = 1 - 2sin(7r/18)) E3 and Archimedean (3, 12 2 ) lat- 
tice (p c = ■JX — 2 sin(7r/18)) EH — the main stream of 
investigation is computational (for analytical approach 
please refer to 17]). For example, for square lattice with 
the nearest neighbours interaction the numerical estima- 
tion of p c (2N) is e = 0.5927460(5) El- 

In this paper the square lattice random site percolation 
is reconsidered. With computer simulation we estimate 
the percolation thresholds for complex neighbourhoods 
Af far beyond standard von Neumann's and Moore's ones. 
We systematically examine neighbourhoods with radius 
r = 1, a/2, 2, 2V2, i.e. which include the 4-th (5N, 
r = 2^2) and the 5-th (6N, r = \/5) nearest neighbours 
of sites. These two neighbourhoods are combined with 
all neighbourhoods with smaller radius r, i.e. the next- 
next-nearest neighbours (4N, r = 2), the next-nearest 
neighbours (3N, r = y/2) and the nearest neighbours (2N, 
r = 1). Some of these percolation thresholds were calcu- 
lated earlier in Refs. [19, 2(jl . As in literature known to 
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us also percolation threshold for 4N+3N neighbourhoods 
is missing we additionally compute its value as well. The 
mentioned above neighbourhoods are presented in Fig. 

ID 



II. CALCULATIONS 

We use Hoshen-Kopelman algorithm [1^ for the oc- 
cupied sites labelling. In the Hoshen-Kopelman scheme 
each site has one label: all sites in given cluster have the 
same labels and different clusters have assigned different 
labels. The labelled lattice examples for 5N+2N neigh- 
bourhoods for the occupation probability p below (a) and 
above (b) the percolation threshold p c are presented in 
Fig. 1 

The sigmoidal curve P{p) tends to the Heaviside func- 
tion H(p — p c ) for L — > 00. The percolation threshold p c 
is evaluated numerically for a finite system as a crossing- 
point of the percolation probability P dependence on 
sites occupation probability p for various system sizes 
L [21( . Such strategy was applied successfully to investi- 
gat e phase transition phenomena, including percolation 
[20l |22| or opinion dynamics [23| . 

Here, basing on P(p) dependence for L = 100, 500 and 
1000 we look for the interval of the length Ap = 10~ 3 
where these three curves cross each other. For example 
for 2N case this interval is (0.592, 0.593). These three sig- 
nificant figures of the left border of this interval estimate 
percolation threshold p c with three digits accuracy. 

In Fig. [3] examples of P(p) dependence for 5N+2N 
neighbourhoods and for L — 100, 500 and 1000 are pre- 
sented. 



III. RESULTS 

The calculated percolation thresholds are collected in 
Tabs. U and [TTJ In Tab. [T] neighbourhoods are sorted 
according to increasing coordination number z while in 
Tab. |TT] according to neighbourhood radius r (and lex- 
icographically). For completeness, in both tables, the 
results from Ref. [2(| are included. The results are also 
presented in Fig. 0] 
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FIG. 1: All square lattice neighbourhoods Af with 5N and 6N: 5N, 5N+2N, 5N+3N+2N, 5N+4N+2N, 5N+4N+3N+2N, 5N+3N, 
5N+4N+3N, 5N+4N, 6N, 6N+2N, 6N+4N, 6N+4N+2N, 6N+5N+2N, 6N+5N+4N+2N, 6N+3N+2N, 6N+3N, 6N+4N+2N, 
6N+5N+2N, 6N+5N+4N+2N, 6N+5N+4N, 6N+5N, 6N+4N+3N+2N, 6N+5N+3N+2N, 6N+5N+4N+3N+2N. Additionally 
the 4N+3N neighbourhoods is presented. With the HKA on a square lattice when we assign the labels for the investigated 
central site (black) we need to check already labelled and occupied sites in its neighbourhood (red slashed sites). The possible 
links to remaining sites in the neighbourhood (black backslashed sites) may be checked later, basing on the neighbourhood's 
point symmetry. 



IV. CONCLUSIONS 

Concluding, we have presented the random site per- 
colation thresholds p c for square lattice with complex 
neighbourhoods which is constituted by 4-th and 5-th 



nearest neighbours. Single site in such neighbourhoods 
has from z — 4 to z — 24 neighbours at distances from 
r = 1 to r = 2V2 from the neighbourhood centre. For 
fixed values of either z or r the percolation threshold 
form a wide interval of values. These ranges overlap each 
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TABLE I: Percolation threshold p c for square lattice site 
percolation and various neighbourhoods constructed with 5N 
an 6N. 



z 


\ f 

vv 


Pc 
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OTVT OM A~\T £?T\T 

2JN , 31N , 4JN , DlN 


/) 

a 


8 


OAT i n\T A AT i OAT /i\T i ,1 AT 

3N+2N, 4N+3N, 6N+4N 


1 — 8 




a AT i OAT /"'AT i OAT 

4N+2N, dJN+oJN 


0.337 • ■ • 






0.270 • ■ • 




/>AT i OAT 

6N+2N 


0.277 ■ ■ • 


12 


/(AT i OAT i OAT /1AT i ^AT i OAT 

4N+3N+2N, 6N+4N+3N 


0.288 • • • 




CAT i OAT 

5N+2N 


0.236 • ■ ■ 




KAT i OAT 

5N+3N 


0.225 • ■ • 




r" AT i A AT 

5N+4N 


0.221 • ■ • 
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6N+3N+2N 


0.240 • ■ • 
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6N+4N+2N 


0.233 • ■ • 




i">AT i CAT 

6N+5N 


0.199 ■ ■ ■ 


16 


CAT i OAT i OAT 

5N+3N+2N 


0.219 • • • 




CAT < ^AT i OAT 

5N+4N+2N 


0.208 • ■ • 




r AT , JAT i OAT 

5N+4N+3N 


0.202 • ■ ■ 




i">AT i CAT i OAT 

6N+5N+2N 


0.187 • ■ • 




6N+5N+3N 


0.182- ■ ■ 




6N+5N+4N 


0.179- ■ ■ 




6N+4N+3N+2N 


0.208- ■• 


20 


5N+4N+3N+2N 


0.196- - ■ 




6N+5N+3N+2N 


0.177- • • 




6N+5N+4N+2N 


0.172- - ■ 




6N+5N+4N+3N 


0.167- ■■ 


24 


6N+5N+4N+3N+2N 


0.164- ■• 



FIG. 2: Example of labelled lattice with 5N+2N neighbour- 
hood for (a) p = 0.20 < p c and (b) p = 0.25 > p c - Three stars 
*** denote one empty site. 



5N+2N: N nm =1000 




0.22 0.225 0.23 0.235 0.24 
occupation probability p 



0.245 



0.25 



FIG. 3: Percolation probability P vs. site occupation prob- 
ability p for 5N+2N neighbourhoods and three lattice sizes 
L — 100, 500 and 1000. The results are averaged out over 
N Iun = 1000 lattice realisations. The common cross-point 
predicts percolation threshold p c . 



other supply earlier observations that p c have not to be 
decreasing function of coordination number z (25j. We 
show the same for neighbourhood radius r. 

As it was stated in Ref. [H p c (6N) = p c (4N) = 
p c (3N) = p c (2N) = 8. Geometrical arguments similar 
to those presented in Ref. [2(| yield 

p c (6N+4N) = p c (4N+3N) = p c (3N+2N), (la) 



and 



p c (4N+2N) = p c (6N+3N) (lb) 



p c (6N+4N+3N) =p c (4N+3N+2N). (lc) 



We improve value of p c (5N) ~ 0.270, as well. 

Our results may support further searching of univer- 
sal formula for p c dependence on space dimension d, site 
coordination number z and additionally the neighbour- 
hoods radius r. However — as we can see from Tab. [TT] 
— these three numbers may be still insufficient to build 
the universal formula prl uK . l2c| as several lattices with 
different neighbourhoods has exactly the same d, z and 
r but different p c . 
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FIG. 4: Percolation threshold p c for various (a) neighbour- 
hood radius r and (b) coordination number z. 
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